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Abstract. In recent years the investigation of hadron structure using lattice tech- 
niques has attracted growing attention. The computation of several important 
quantities has become feasible. Furthermore, theoretical developments as well as 
progress in algorithms and an increase in computing resources have contributed 
to a significantly improved control of systematic errors. In this article we give an 
overview on the work that has been carried out in the framework of the Hadron 
Physics 13 (I3HP) network "Computational (lattice) hadron physics". Here we 
will not restrict ourselves to spin physics but focus on results for nucleon spec- 
trum and structure from the QCDSF collaboration. For a broader overview of 
developments in this field see, e.g., [l]. 

1 Introduction 

The main goal of lattice hadron phenomenology is the computation of generalized parton distri- 
butions (GPDs) of mesons and nucleons from first principles. GPDs have become a theoretical 
framework which enables us to study many fundamental aspects of the intrinsic hadron struc- 
ture. It allows us to confront experimental results and lattice calculations since it includes 
as limiting cases hadron form factors as well as polarized and unpolarized parton densities, 
i.e. quantities which are investigated by various experiments at CERN (COMPASS), DESY 
(HI, Zeus, Hermes) and JLab. Combining information from both experiment and lattice offers 
new opportunities to explore hadrons as extended objects. 

On the lattice we are restricted to the calculation of moments of GPDs. As an example we 
consider the unpolarized GPDs H'^{x,£,,Q^) and E'^{x,£,,Q'^). Their moments are related to 
the generalised form factors (GFFs) Al ^{Q^), Bl ^{Q'^) and C^{Q'^): 
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Fig. 1. The left plot shows a comparison of Zy obtained from 1-loop perturbation theory (dashed 
line) and a non-perturbative determination (solid line). The right plot shows the values for ro(/3,mps) 
(/3 = 5.20, 5.25, 5.29, 5.40) together with an extrapolation to the chiral limit. 



The GFFs can be calculated from a parametrisation of (nuclcon) matrix elements. As an illus- 
trative example one may consider the GFFs A^ q and B\q: 



1v 



(p',s'|u7^^|p,s) = ii(p',s')^7^<onQ')+«'^'^''^sS(Q') \u{v.s), 



(2) 



where p (s) and p' (s') denote initial and final momenta (spins), q = p' — p the momentum 
transfer (with Q^ — —q^). Instead of calculating the matrix elements (i.e., the l.h.s. of Eq. ^) 
in Minkowski space-time, we calculate these within a discretized Euclidean space-time by using 
appropriate combinations of two- and three-point correlation functions. A typical combination 
is the following ratio: 
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(3) 



where C2{t,p) is the unpolarized two-point function with a source at time and sink at time 
t, while the three-point function C3(t, r,p',p) has an (local) operator O insertion at time r. 
The ratio R{t,T,p' ,p) is expected to become constant when keeping sufficiently far away from 
source and sink, i.e. for <C t ^ i < 5-^r- For further details see, e.g., [2]. 



1.1 Lattice obstacles 



All lattice calculations are faced with the problem of keeping their systematic errors under 
control. The main sources for such systematic errors are the following: 

— In almost all cases lattice results have to be extrapolated to the 

— infinite volume limit: 1^ ^ oo. 

— continuum limit: a — > where a is the lattice spacing. 

— region of light quark masses: nips ~^ it^tv- 

— Renormalization and mixing of operators. 

— Conversion of the lattice results into physical units. 

The results which we present in this article have been obtained from simulations with Nf — 2 
flavours of dynamical 0(a)-improved Wilson fermions (so-called Clover fermions). Improvement 
of both the action and the operators is expected to bring discretization errors down to 0(0^). 
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Fig. 2. The left plot shows lattice results for ttin together with a fit to Eq. Q. The star in- 
dicates the physical point. In the right plot we compare the relative finite size effect S{L) = 
{mN{L) — mjv(oo))/r7i]v(oo) as measured on the lattice and predicted by BxPT. The diamonds show 
results from lattices with F = 12^ x 32, 16^ x 32, 24^ x 48 (left to right). 



Our simulations have been performed with lattice spacings in the range of 0.11 down to 0.07 fm. 
For most of our results we find discretization effects to be small and often negligible compared 
to, e.g., statistical errors. 

For some of the observables considered here, finite volume effects turned out to be large. In 
our simulations the lattices have a spatial extension in the range L^ = 1.4, ■ • ■ , 2.6 fm. Results 
from chiral effective field theories may be used to estimate or even correct for finite size effects. 

Typically, the largest uncertainty originates from the extrapolation of the results down to 
the region of physical light quark masses. Results presented here have been obtained for pseudo- 
scalar meson masses in the range of 300 MeV up to 1 GeV. Due to a combination of significantly 
improved algorithms (see [3] for a recent review on algorithms for simulations with dynamical 
fermions) and the availability of increased computer power it has become possible to start 
exploring a quark mass region where one may hope results from chiral perturbation theory to 
be applicable. 

Another source of systematic errors stems from renormalization. Perturbative renormaliza- 
tion techniques suffer from their complexity when going beyond 1-loop and, more importantly, 
from their poor convergence. It turned out to be crucial to employ non-perturbative renormaliza- 
tion techniques. As an illustrative example we compare in Fig. [1] (left panel) the vector-current 
renormalization constant Zy obtained from 1-loop perturbation theory and a non-perturbative 
determination, where we imposed as renormalization condition the equivalence of the local and 
the conserved vector current (see U for a detailed description of this method) . 

Let us finally consider the problem of converting lattice results into physical units. It has 
become popular to use the Sommer parameter vq, which is defined by the expression Tq F{ro) — 
1.65. The force F{r) can be determined with rather high statistical accuracy from the static 
quark potential. This strategy however suffers from the problem that the phenomenological 
value is not very well known. Typically, ro ~ 0.5 fm is used. Lattice results, e.g., for the nucleon 
mass (see [5]), suggest a significantly smaller value. Here we will use tq — 0.467 fm [617] . Fig. [1] 
(right panel) shows our results for tq. 



2 Masses 



As a first example for using chiral effective field theories in order to guide extrapolation to 
physical quark masses we consider the mass of the nucleon, ttin ■ Calculations using relativistic 
baryon chiral perturbation theory (B^PT) suggest a rather non-trivial quark mass dependence. 
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Fig. 3. The left plot shows ^^(m-ps) in the infinite volume limit and the right plot the data points 
before correcting the finite size effects. 

The p-expansion at 0(p*) in an infinite volume gives the following quark mass dependence [8]: 
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Even using all our lattice measurements at quark masses in the the range < Trips ^ 
650 MeV it is not possible to sufficiently constrain all parameters in Eq. (Uj). We therefore 
reduce the number of free fit parameters to the nucleon mass in the chiral limit Mq, the not 
very well known low-energy constant (LEG) ci and the counter-term ei{X) (we use A = 1 GeV). 
For all other parameters, i.e. the LEGs C2 and C3, the pion decay constant Fq and the nucleon 
axial coupling gA,o, we use phenomenological estimates. Both lattice data and the resulting fit 
are shown in Fig. [2] (left panel). The lattice data seems to fall on a universal curve indicating 
discretisation effects to be small, which we thus ignored. We observe that m^(mps = TOtt) 
is consistent with experiment. Furthermore, we find ci = —1.02(7) GeV~^, a value which is 
consistent with other estimates [9]. 

BxPT has also been used to calculate the finite size effects 0: 



mN{L,) - TOjv(oo) = Aa{L,) + Ai,{L,) + 0{p^), 



(5) 



where Aa represents the 0{p^) result for the volume dependence and Z\b additional contribu- 
tions at 0{p^). All coefficients in Aa{Ls) and Ai,{Ls) are also present in Eq. (j4|). It is thus 
possible to estimate the finite volume effects obtained from Eq. ([5]) using the coefficients ob- 
tained from a fit to Eq. ((4]). These estimates can then be used for a comparison with values 
of the nucleon mass computed on lattices that have a smaller (physical) volume. Fig. [2] (right 
panel) shows such a comparison at fixed m-ps — 590 MeV. 



3 Nucleon axial coupling 

Next, we consider the axial couphng constant qa = G'a(O), where Ga{Q^) is the axial form 
factor of the nucleon. It is determined from the renormalized axial vector current A^ = Za (1 + 
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bA amq)A^, where amq is the bare quark mass. Here we only look into the iso- vector case where 
contributions of so-called disconnected terms cancel. These are contributions of diagrams which 
only interact with the hadron via the exchange of gluons and which are hard to compute on 
the lattice. While Za is known non-perturbatively jlOj , bA is only known perturbatively and is 
computed using tadpole improved one-loop perturbation theory. 

Like for the nucleon mass the quark mass dependence of the iso- vector nucleon axial cou- 
pling has been calculated using chiral effective field theory |ll|10j . The calculations using the 
small scale expansion (SSE) have been performed in the infinite volume limit as well as for a 
finite spatial cubic box of length Lg- We define gAirrips, L^) — gA{mps) + A9A(m-ps, Lg), where 
AgA{mpsT Ls) denotes the finite size effects. Given the large number of parameters in the re- 
sulting expressions we again have to fix some of them using phenomenological input like the 
pion decay constant, the NA mass splitting and the axial NA coupling. 

We find the finite size effects to be significantly larger compared to the nucleon mass. These 
effects can not be ignored and we have to resort to a fit which includes both quark mass 
dependence and finite size effects. In Fig. [3] we show results for the quark mass dependence in 
the infinite volume limit and finite size effects. 

For results for Ga{Q^) at non-zero Q^ and for the pseudo-scalar form factor Gp{Q^) see 
[l2] . For an independent calculation of qa see ll3j . 



4 Electro-magnetic form factors 

JLab polarisation experiments [14I15J have in recent years led to a revived interest in the nucleon 
electro-magnetic form factors. Measurements of the ratio of the proton electric to magnetic form 

factors, z^^P-'Ge (Q^)/Gm (Q^), showed an unexpected decrease. This means that the proton's 
electric form factor falls off faster than the magnetic form factor. 

The form factors are obtained from the standard decomposition of the nucleon electromag- 
netic matrix elements 



{p\s'\J^\p,s) = u{p,s) 



Zmpf 



u{p, s) . (6) 



By calculating the matrix elements on the l.h.s. and the nucleon mass we obtain the Dirac form 
factor Fi{Q'^) and the PauH form factor i^2(Q^)- 

Particularly interesting are the Q^ scaling and flavour dependence of the form factors. We 
find that the lattice data can be well parametrized using a pole ansatz 

^^^^'^ = [i + Q^/Mfy ^^) 

(For a different ansatz see [T2].) Naively, one would expect p = 2 for Fi and p — i for i^2- 
Our lattice data, however, has a flavour dependence which favours p = 2 for i^^ and p = i 

otherwise. In Fig. 0] (left panel) we plot the ratio F{ / F{ which one naively would expect to 
be constant. This result should be taken with some care since disconnected contributions have 
not been calculated. These contributions only cancel in the iso- vector case. Our observation is 
however consistent with the flavour dependence observed in fits to experimental data (TB]. 

For a more quantitative comparison with experimental data we have to extrapolate our 
results to physical quark masses. We again do this on the basis of results from chiral effective 
field theories using the SSE [17118] . These calculations predict a strong quark mass dependence 
for the iso- vector form factor radii and the iso- vector anomalous magnetic moment in the small 
quark mass region. In Fig. |4] (right panel) we compare our results for the iso-vector Dirac 

radius Tj and a fit to the SSE expression. The lattice results are signiflcantly smaller than 

the experimental value. However, the SSE results indicate a strong quark mass dependence for 
TOps < 250 MeV, a region which is currently difficult to access in lattice simulations. 
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4. The left plot shows the ratio F^ /F{' as a function of Q for different /3 values but similar 
400 MeV. The right plot shows the extrapolation of our results for the Dirac form factor radius 
The solid line shows a fit with the SSE expression, the dashed line a fit to an ansatz linear in 



the quark mass. 



5 Moments of unpolarized structure functions 

As another quantity which can be calculated on the lattice we consider the lowest moment of 
the unpolarized nucleon structure function, (x) = ^20(0), where Ag q is the first moment of 

This moment is determined from the matrix element 



(P, . 



q{i^iP-2} _ fpj. 



\p,s)^2A%[p^^p'^--Tt] 



(8) 



where O^^'^^^' = wy^'^iD^'^u. The renormalisation of this matrix element has been performed 
non-perturbatively by means of the RF-MOM method [19]. In the perturbative conversion to 

MS at a scale of 2 GeV we have used A^^ = 261 (17) (26) MeV [20]. 

In the range of quark masses which is currently accessible the results in the iso- vector channel 
^,j,^(u-d) g^j.g significantly larger than the experimental value. However, it has been suggested that 
this quantity may become significantly smaller at very light quark masses j21| . This has been 
confirmed by recent calculations in the framework of baryon chiral perturbation theory (B^PT) 
[22j . These recent calculations have been used to fit our results using the pion decay constant, 
the nucleon mass and the nucleon axial coupling in the chiral limit as phenomenological input. 
This leaves two free fit parameters, which allows us to fit the data for different values of /3 
separately restricting the fit range to < irips < 650 MeV. In Fig. [5] we plot the fit to the 
(3 = 5.29 data showing little evidence for (a;) '""''-' (mps) becoming smaller for nips -^ rn-,^. In our 
results we find no indication for large discretisation effects (see right panel of Fig. [5]). However, 
for the results at light quark masses finite size effects cannot be excluded. 

In Fig. [5] (left) the Q^ dependence of the generalized form factor A\ q is shown together with 
the other unpolarized GFFs Bf ^-i^d Cf . The lattice results for A2 g can be well described by a 
dipol expression A^QiQ"^) — ^2.o(0)/(l + Q^ A^l))^- The quark mass dependence of the dipole 
mass is also shown in Fig. [5] (right plot). For further details see [23] . 



6 Other results 



Several other results have been published in recent years, which we will summarize here only 
very briefly: 
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Fig. 5. The left plot shows our results for (s)*^" '*'(?nps) at /3 = 5.29 (the results for other val- 
ues of /3 are similar). The solid line denotes a fit to a BxPT expression [22]. The right plot shows 
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as a function of the lattice spacing together with a fit to a constant. The star 



denotes the phenomenological value from CTEQ6. 
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Fig. 6. The left plot shows the generalised form factors Aj g, B|o and C2 in the iso- vector channel at 
K — 0.13632 and /3 = 5.29 with rrips ~ 350 MeV. In the right plot the dipole masses of A2 q are shown 
as a function of the pseudo-scalar meson mass. The star represents the experimental value of the tensor 
meson mass /2. 



A calculation of the lowest two moments of transverse spin densities of quarks in the nucleon 
has revealed strongly distorted densities of transversely polarized quarks in the nucleon 24J . 
The electro-magnetic form factor F^ of the pion calculated on the lattice has been found 
to be in good agreement with experimental results |25j . The charge radius squared was 
obtained as (r^) = 0.441(19) fm^ 

Lattice calculations of the first two moments of the quark tensor GPD EJ^ indicate - like in 
the case of the nucleon - a strongly distorted spatial distribution of the quarks if they are 
transversely polarized |26j . 



7 Conclusions and outlook 



In this talk we have presented a small selection of results to demonstrate where lattice calcula- 
tions can contribute to the investigation of hadron structure. While these results may improve 
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our qualitative understanding, for quantitatively precise results further efforts are needed and 
keeping systematic errors under control remains a major challenge. In particular, the necessary 
extrapolations to the infinite volume, vanishing lattice spacing and to the physical quark masses 
remain a major source of uncertainties. While results from chiral effective field theories help 
to guide these extrapolations, eventually simulations using large lattices and very light quark 
masses are needed. We conclude with the good news that such simulations are now becoming 
feasible. 
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